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Abstract. We present near-to-mid-infrared spectral energy distributions (SEDs) for 21 Seyfert 
galaxies, using subarcsecond resolution imaging data. Our aim is to compare the properties 
Seyfert 1 (Syl) and Seyfert 2 (Sy2) tori using clumpy torus models and a Bayesian approach 
to fit the infrared (IR) nuclear SEDs. These dusty tori have physical sizes smaller than 6 pc 
radius, as derived from our fits. Active galactic nuclei (AGN) unification schemes account for 
a variety of observational differences in terms of viewing geometry. However, we find evidence 
that strong unification may not hold, and that the immediate dusty surroundings of Syl and 
Sy2 nuclei are intrinsically different. The Type 2 tori studied here are broader, have more 
clumps, and these clumps have lower optical depths than those of Type 1 tori. The larger the 
covering factor of the torus, the smaller the probability of having direct view of the AGN, and 
vice-versa. In our sample, Sy2 tori have larger covering factors (Ct=0.95±0.02) and smaller 
escape probabilities than those of Syl (Ct=0.5±0.1). Thus, on the basis of the results presented 
here, the classification of a Seyfert galaxy may depend more on the intrinsic properties of the 
torus rather than on its mere inclination, in contradiction with the simplest unification model. 



1. Introduction 

The unified model for active galaxies [I] is based on the existence of a dusty toroidal structure 
surrounding the central region of AGN. This geometry of the obscuring material allows the 
central engines of Type 1 AGN to be seen directly, resulting in typical spectra with both narrow 
and broad emission lines, whereas in Type 2 AGN the broad-line region (BLR) is obscured. 

Pioneering work in modelling the dusty torus [281 123 EH HI E] assumed a uniform dust density 
distribution, although from the start, [16J realized that smooth dust distributions cannot survive 
within the AGN vicinity. In order to solve the discrepancies between the observations and the 
predictions from smooth torus models, an intensive search for an alternative torus geometry has 
been carried out in the last decade. The clumpy dusty torus models [22], [231 021 G2] propose 
that the dust is distributed in clumps, instead of homogeneously filling the torus volume, and 
they are making significant progress in accounting for the mid-IR (MIR) emission of AGNs 

[E1EIJE1E51I331EDII21II]. 



In this work, we report IR nuclear fluxes for a representative sample of nearby Seyfert galaxies, 
for which we estimate unresolved nuclear MIR fluxes. We compile near-IR (NIR) nuclear fluxes 
from the literature of similar resolution to construct nuclear SEDs for all the galaxies, and 
fit them with clumpy torus models, which we interpolate from the |22[ [23] database. Table 1 
summarizes key observational properties of the sources in the sample. 

Table 1. Basic Galaxy Data. 
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2. MIR Observations 

Ground-based MIR high-angular resolution observations of 21 nearby active galaxies were carried 
out over the past years for a variety of science drivers. We make use of this archive of data in 
this paper. The first set of observations was obtained with the University of Florida MIR 
camera/spectrometer OSCIR at the 4 m Blanco Telescope at Cerro Tololo Inter- American 
Observatory (CTIO) and at the Gemini North Telescope. Another set of observations was 
performed with the MIR camera/spectrograph T-ReCS (Thermal- Region Camera Spectrograph) 
on the Gemini-South telescope. The last set was obtained with the MIR camera/spectrograph 
Michelle on the Gemini North Telescope (see Table 1). Observations of PSF and flux standard 
stars were made for each galaxy through the same filters. The PSF star images were employed 
to determine the unresolved (i.e., nuclear) component of each galaxy. For further details on the 
MIR observations of the whole sample we refer the reader to |33[ [31] . 

3. SED Observational Properties 

Using the MIR data described in Section [2] and NIR data of similar resolution compiled from the 
literature, we constructed subarcsecond resolution nuclear SEDs in the wavelength range from 



~1 to 18 /jm for the galaxies in Table 1. Figures [T] and [2] show the individual SEDs. 

Using the highest resolution data only (<0.55 arcsec) we constructed average Syl and Sy2 
templates (see Figure 4 in [31]) and measured the 1.265-18.3 fim IR slope (f„ a i/~ aiR ) of the 
Syl template: am = 1.7 ±0.3, as well as the NIR (a^m = 1.6 ±0.2) and MIR spectral indexes 
{oimir = 2.0 ± 0.2). A flat NIR slope indicates an important contribution of hot dust emission 
(up to ~1000-1200 K; [6J) that comes from the immediate vicinity of the AGN. The shape of the 
Sy2 mean SED is very steep {am = 3.1±0.9, a^iR = 3.6±0.8, and aMiR = 2.0±0.2) compared 
with those of the Type 1 Seyferts. In general, Sy2 have steeper 1-10 fim SEDs than Syl [HI El [2]. 
On the contrary, the MIR slope is the same for the Syl and Sy2 templates [oimir = 2.0 ± 0.2). 

4. SED Modelling 

The clumpy dusty torus models of Nenkova et al. (2002) hold that the dust surrounding the 
central engine of an AGN is distributed in clumps, instead of homogeneously filling the torus 
volume. The clumpy database now contains 1.2xl0 6 models, calculated for a fine grid of model 
parameters. The inherent degeneracy between these parameters has to be taken into account 
when fitting the observables. To this end, we developed a Bayesian inference tool (BayesClumpy; 

0)- 

The priors for the model parameters are assumed to be truncated uniform distributions in 
the intervals reported in Table 2. Therefore, we give the same weight to all the values in each 
interval. Apart from the six parameters that characterize the models, there is an additional 
parameter that accounts for the vertical displacement required to match the fluxes of a chosen 
model to an observed SED. This vertical shift scales with the AGN bolometric luminosity (see 
Section [5j). 



Table 2. Clumpy Model Parameters and Considered Intervals 



Parameter 


Abbreviation 
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Width of the angular distribution of clouds 
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Radial extent of the torus 


Y 
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Number of clouds along the radial equatorial direction 
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[1, 15] 


Power-law index of the radial density profile 
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Inclination angle of the torus 
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Optical depth per single cloud 
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The results of the fitting process of the IR SEDs with the interpolated version of the clumpy 
models of [2'2\ are the posterior distributions for the six free parameters that describe the 
models and the vertical shift. When the observed data introduce sufficient information into 
the fit, the resulting posteriors will clearly differ from the input uniform priors, either showing 
trends or being centered at certain values within the intervals considered. 

We fitted the individual IR SEDs with BayesClumpy, modelling the torus emission and 
the direct AGN contribution (the latter as a broken power law). We also consider the IR 
extinction curve of [7] to take into account any possible foreground extinction from the host 
galaxy. Although the solutions to the Bayesian inference problem are the posterior distributions 
of each parameter, we can translate the results into corresponding spectra (see Figures Q] and 
[2]). The solid lines correspond to the model described by the combination of parameters that 
maximizes their probability distributions (maximum-a-posteriori; MAP). Dashed lines represent 
the model computed with the median value of the probability distribution of each parameter. 



Shaded regions indicate the range of models compatible with the 68% confidence interval for 
each parameter around the median. 
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Figure 1. High spatial resolution IR SEDs of the Seyfert galaxies (dots). Solid and dashed lines correspond 
to the MAP and median models, respectively. Shaded regions indicate the range of models compatible with the 
68% confidence interval for each parameter around the median. 

The clumpy models successfully reproduce the observed Seyfert SEDs studied here with 
compatible results among them. This indicates that the NIR and MIR unresolved fluxes 
employed here are dominated by a combination of reprocessed emission from dust in the torus 
and direct AGN emission. 
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Figure 2. Same as in Figure [T] but for the remaining Sy2 galaxies. 
5. Comparison between Type 1 and Type 2 Seyfert tori. 

The aim of this work is to compare between the properties of Type 1 and 2 tori under the 
assumption that the SEDs studied here are torus/ AGN dominated. Despite the relatively low 
number of objects considered (7 Syl and 9 Sy2; the unreliable fits of NGC 1808 and NGC 7582 
are excluded), we find that some of the parameters are significantly different between Syl and 
Sy2. 

To take full advantage of the Bayesian approach, the best way to compare the results for Syl 
and Sy2 galaxies is to derive joint posterior distributions for the full Type 1 and Type 2 datasets 
respectively. Thus, we normalized all the Syl SEDs at 8.74 /im and fitted them together using 
BayesClumpy, and we did the same for the Sy2. We considered the mean redshift for the Syl 
(z=0.0061±0.0045) and for the Sy2 (0.0078±0.0051) in the fits. In Figure [3] we show the Syl 
(left panel) and Sy2 fits (right panel). Note that the MAP and median models predict a flat 
SED with the silicate feature in weak emission for the Syl galaxies, and steeper and with the 
silicate band in shallow absorption for Sy2. 

The comparison between the Syl and Sy2 posterior distributions is shown in Figured! From 
a visual inspection it is clear that the joint posteriors of the parameters No, q, ry, and o are 
completely different between Syl and Sy2. There is no overlap between the 1-sigma intervals. 
In Table 3 we report the median and mode values of the histograms in Figure H] 

In order to quantify how different the probability distributions are, we calculated the 
Kullback-Leibler divergence (KLD; |17j ) between the Syl and Sy2 posteriors. This divergence 
takes into account the full shape of the posterior and it is a positive value, and equal to zero 
when two distributions are identical. Therefore, the larger the value of KLD, the more different 
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Figure 3. Same as in Figure Q] but for the Syl (left) and Sy2 SEDs (right) normalized at 8.74 pm. 




Figure 4. Posterior distributions resulting from the fits of the joint Syl and Sy2 SEDs. KLD values derived 
from the comparison between Syl and Sy2 for each parameter are labelled. 



Table 3. Statistics of the comparison between Syl and Sy2 parameters 
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the posteriors. We find KLD>1 for a, No, q, and ry (see Figure H]). These are indeed the four 
parameters whose 1-sigma regions do not overlap. Thus, we consider their differences significant 
between Syl and Sy2. For both Y and i we find KLD<1 and similar median values between 
Syl and Sy2. 

Syl tori are narrower and have fewer clouds (a = 44°±7 ; No=4±l) than those of Sy2 
(a = 63°±| ; No=ll±^). The radial density distribution of the clouds is also different between 
the two Seyfert types according to this analysis: in Sy2, the majority of the clumps are 
distributed very close to the nucleus (i.e. steep radial density distribution; g=2.3±0.1) whereas 
for Syl the clouds distribution is flatter (g=0.8±0.2). On the other hand, the optical depth of 
the clouds in Syl tori is larger (ry=133±g) than in Sy2 (ry=30±l). 



Interestingly, we find high as well as low values of the inclination angle of the torus for Syl 
and Sy2 (see Table 3). This variety in the i values translates into the similar median values of 
the joint Syl and Sy2 posteriors (47°±go for Syl and 54°±i5° f° r Sy2), which are intermediate 
within the considered prior (i=[0°,90°]). This is telling us that, in the clumpy torus scenario, 
the classification of a Seyfert galaxy as a Type 1 or 2 may depend on the intrinsic properties of 
the torus rather than in mere inclination. 

In Figure \5\ we represent the median values of a and No for the different Seyfert types over 
the covering factor contours. The covering factor is defined as Ct = 1 — J e~ NLOS ^dcos(i). 
Type 1 nuclei tend to be located within lower Ct contours (Ct <0.6) than those of Type 2s, for 
which Ct >0.5, with the exceptions of Centaurus A and Mrk 573. We have represented with 
larger symbols the median values from the joint a and Nq posteriors reported in Table 3. 

Since the covering factor is a non-linear function of the torus parameters, we took full 
advantage of our Bayesian approach and generated joint posterior distributions for Ct from those 
in Figured The median values of the histograms are CT(Sy2)=0.95±0.02 and Cr(Syl)=0.5±0.1. 
The divergence between the Syl and Sy2 Ct posteriors is KLD=28, indicating that the difference 
is significant (1-sigma regions do not overlap). Thus, Syl tori in our sample have lower Cts 
than those of Sy2, implying that they are intrinsically different. 
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Figure 5. a versus No for the individual galaxies. Either median values or upper/lower limits are taken from 
the fits presented here. Error bars indicate 68% confidence level around the median. The big dot and square 
correspond to the average a and No values for Syl and Sy2 from Table 3. 

The clumpy model fits yield the intrinsic bolometric luminosity of AGN (L^ N ) by means of 
the vertical shift applied to match the observational data points. Combining this value with the 
torus luminosity (L^ ), obtained by integrating the corresponding model torus emission (without 
the AGN contribution), we derive the reprocessing efficiency (RE) of the torus (V^/L^ N ). 
The previous values are calculated on the Bayesian framework, by combining the posterior 
distributions of the model parameters. Sy2 tori in our sample are more efficient reprocessors 
than Syl, absorbing and re-emitting the majority of the intrinsic AGN luminosity in the IR: 
RE(Sy2)=[0.4, 1.0], with a median value of 0.8 and RE(Syl)=[0.2,0.7], with median of 0.5. 

We considered a possible dependency of the RE (or alternatively the covering factor; Ct) on 
Lf^f N , since the amount of incoming radiation from the AGN could possibly have some influence 
on the reprocessed energy or even in the torus properties (e.g., receeding torus scenario; |18j). 
However, we find no relationship between the two quantities in the luminosity range considered. 



This means that the reprocessing efficiency depends primarily on the total number of clouds 
available to absorb the incident radiation, i.e. on the torus covering factor. However, the possible 
dependence of the torus properties on the AGN luminosity considering a broader luminosity 
range is further investigated in JTj. 

In general, the IR SED fitting does not constrain the size of the torus (Y) as well as other 
model parameters (see Section 5.3 in [33] )■ The NIR and MIR observations are sensitive to 
the warm dust (located within ~10 pc of the nucleus), which depends on the combination 
of model parameters No, q, and Y [35]. Far-IR observations are more sensitive to the torus 
extent independently, as shown in [32]. The outer size of the torus scales with the AGN 
bolometric luminosity: R a = YR^, so assuming a dust sublimation temperature of 1500 K, 
R a = 0.4 Y (L^^/IO 45 ) 0,5 pc. We derived R Q posterior distributions from those of and 
Y and find that all tori in our sample have outer radii smaller than 6 pc, in agreement with MIR 
direct imaging of nearby Seyferts [26, 30J and also interferometric observations [T5 | |36" 1 [20 ], [29] . 

Summarizing, we find tantalizing evidence, albeit for a small sample of Seyfert galaxies and 
under the clumpy torus hypotesis, that the classification as a Type 1 and Type 2 depends on the 
intrinsic torus properties, in contradiction with the simplest unification model. 
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